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1. Introduction 



The Yang-Baxter (YB) equations assure integrability in two-dimensional lattice 
models and quantum field theories. The YB equations have the form 

R 12 (9' - 9)R 13 (9)R 23 (9') = R 2 3(9')R 13 (9)R 12 (9' - 9) , (1.1) 

where the i?-matrix elements R^(9) (with 1 < a,b,c,d < q , q > 2) define the 
statistical (or Boltzman) weights for a vertex model in two dimensions. The general 
Non Intersecting String (NIS) model is a q(2q — 1) vertex model introduced in [1] 
as a solution of the Yang-Baxter equations. It generalizes the q = 2 ferroelectric 
models and the q = 3 model of Stroganov [2]. 

The non-zero Boltzmann weights of the model are depicted in Fig.l. We distin- 
guish between indices of two different kinds that we shall call "red" and "black" , 
corresponding to dashed and solid lines respectively in the figure. The red in- 
dices, denoted by a,f3, . . . , take on r values, whereas the black ones, denoted by 
fi, u, . . . , take on the remaining t values, with r + t = q > 2 . Latin indices a,b, . . . 
run over all the possible values 1, . . . , q . The values of the R -matrix elements in 
Fig.l are: 

R%%{9) = sinhtf + sinh(?? - 6) , R^{9) = smh V , 

Rf p {9) = sinh^ - 9) , a ± (3 , R%{6) = e e sinh^ - 61) , ^ v , 

RaW = e 6/2 smh(7? - 9) , Rf d (9) = Rf a ( V - 9) . 

(1.2) 

The Boltzmann weights of the NIS model, like the six-vertex ones, depend on the 
spectral parameter 9 and on the anisotropy parameter 77 but, differently from the 
six-vertex case, r\ is not continuous but takes on discrete values, solutions of the 
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second order equation [1] 



l)e 2r > + re* - 1 = . 



1.3) 



We have then the two solutions rj and rf + 27r given by: 



^ _ ^/r 2 + 4(g - r - 1) - r 
~ 2(q-r- 1) 



(1.4) 



with the limiting case e v = 1/r for g = r + 1 and 



+ 4(g _ r _ x) + r 



2(9 



1) 



1.5) 



We see that 7/ > for g < 2(r + 1) and rj' < for q > 2(r + 1) . It must be 
noticed that although the R -matrix becomes complex for the solution rj' + in , 
the imaginary factors cancel out in the partition function for any finite size. For 
simplicity, in the following we shall consider the solution (1.4) which corresponds 
to real negative 77 and real R -matrix. 

By definition, (see last equation in (1.2)) the model is crossing invariant and i] 
is the crossing parameter. As it appears clear from Fig.l, the NIS model describes 
configurations of non-intersecting patterns of polygons on the square lattice. For an 
iV -sites system, the spin Hamiltonian associated to the NIS model can be written 

as 



N 



n(a, b) 



n(a, b) 



o,6=l V / a,6=l 

<? 

- sinh r] n ( a > a ) e T ® e j+l 
0=1 



sinh?] - coshr^^ e°j a <8> ef +1 



where n(a, b) is the number of black indices among a and b ( n(a, b) = 0, 1 or 2 ) 
and {e ab ) cd = 5 ac 5 bd . 
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Several papers dealing with the special case t — , also called the "separable 
vertex model", already appeared in the literature [1-5]. In particular, r = q = 3 
corresponds to the model considered by Stroganov [2]. Equivalence between the 
separable NIS model with a q 2 -Potts model has been established in [4]. Analytical 
continuation to — 2 < q < 2 of the separable case corresponds to the scaling limit 
of the 0(n) vector model with — 2 < n < 2 away from the critical point [5]. As 
it is well known, the partition function of integrable vertex models can be used for 
computing knot invariants [6]. The separable NIS model leads then to the Jones 
polynomial [7]. 

In this paper we find the exact solution of the NIS model in its full generality 
(t 7^ 0) by means of analytic Bethe Ansatz (BA). We show that with periodic 
boundary conditions on a A -sites lattice ( A even), the eigenvalues of the transfer 
matrix r(9) of the general NIS model have the form 

A(0, A) = F(9) smh N (r] - 9) e eB + F( V - 9) sinh^ 9 e^ B , 

where B is an integer ( < B < N ), A = (Ai, . . . , A m ) and 

f(9\ = P ™ n sin ( A i - w - Wg) 

{) ^8^-19 + ^/2) ' 

The numbers Xj in F(9) are obtained by solving the BA equations which read, 
for the NIS model, (see Sect. 3): 

N 



( sin(A fc + ^/2) \ e 2i(\ kB -u) = TT 
\sm(X k -ir]/2)J f = j 



sm(A fc - Xi + irf) 

—— — , l<k<m. (1.6) 

sin(Afc - Ai - it]) 



Here w is a real phase. We conjecture (see Secs.3 and 6 for more details and a 
proof in the limit N — > oo ) that e tuJ is in fact a root of unity i.e. 



e iul 



for some integer I . The eigenvectors of t{9) are thus labelled by B , uj and m . 
The antiferromagnetic ground state is obtained for B = , m = N/2 and uo = . 
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The corresponding eigenvalue of t(6) coincides exactly with the antiferromagnetic 
ground state of the six-vertex model with the identification 7 = —r] . Notice that 
the BA equations (1.6) are identical to those of the six-vertex model if 7 = —77, 
except for the phase factor e M x kB-uj) _ 

We recall that for a spin-s model the L.H.S. of the BA equations is given by 

[10] 

/ sin(A fc + isj) \ Ns 
Vsin(A fc - isj)J 

whose s — > 00 limit is exp(— 2i\} t N s ) . Moreover, as shown in ref. [9], integrable 
models with mixed spins can be constructed and solved by BA. This evidence 
suggested us that a deep link exists between the NIS model and a model contain- 
ing both spins 1/2 and infinite spins. We are in fact able to prove that the BA 
equations and the eigenvalues of the NIS model coincide with those of a model 
containing iV spins 1/2 and N s spins s in the limit in which s is infinite. The 
proof of this equivalence is given is Sec. 6 in the limit N, N s — > 00 . As a byproduct 
we find that the phases e lU) are roots of unity. In this context, the discrete param- 
eter B turns out to be the contribution of the 00 -spins to the total z -component 
of the spin (normalized to one) for each BA state. 

It is known that it is possible to associate a family of solutions of the YB algebra 
(1.1) to each Lie algebra G [8]. Trigonometric (hyperbolic) solutions usually pro- 
vide highest weight or cyclic representations of the corresponding quantum group 
Gq . Here q = e* 7 ( e 7 ) where 7 is the anisotropy parameter in the R -matrix. As 
the discussion above indicates, the NIS model should provide a new representation 
of the quantum group SU(2)q with q = e n (notice that \q\ ^ 1). This repre- 
sentation is absent in the "classical" case q — 1 and is probably indecomposable 
containing s = 1/2 and s = 00 . 

The plan of the paper is as follows. In Sec. 2 we summarize the general 
properties of the NIS R -matrix and transfer matrix. In Sec. 3 we construct the 
BA solutions of the NIS model through analytic BA. Sec. 4 contains checks for 
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large values of the spectral parameter 6 of the results obtained by BA. In Sec. 5 
we find the free energy and the excitations in the thermodynamic limit. Finally, 
in Sec. 6, the connection between the NIS model and a spin 1/2-spin oo model is 
derived. 

2. General properties of the NIS model 

For the general NIS model we find, by direct computation, the following prop- 
erties for the R matrix 

a) Regularity at 9 = 

R(0) = sinhr? 11 . (2.1) 

b) Unitarity 

R(0)R(-0) = p(6) 11 , (2.2) 
where p{6) = sinh(r] — 6) sinh(^ + 6) 

c) Quasi periodicity 

R(e + m) = (u®l)R(6)(u- 1 ®h) , (2.3) 

where u and h are the diagonal matrices = u a 5 a i) and h a i, = h a 5 a i) with 

u a = 1 , u il = % , J 1 <a <r 
h a = -1 , hfj, = 1 . [r + 1 < /i < q 

d) Invariances. A Yang-Baxter algebra is said to be invariant under a transfor- 
mation group Q if the R -matrix satisfies [11] 

[g ® g , R(9)] = , WeC, VgeG . 

In the separable case, r — q , from eqs.(1.2) it follows 

^ T = l , 

that is we have invariance under the orthogonal group. 



6 



In the other limiting case in which all the indices are of the black kind ( r = ), 
we find 

(gg T )w = i , V/i 

and 

which, for the linear independence of the factors e^^ -7 ^ yields 

9w9v P = , /i^ ■ 

Thus Q is formed by matrices g with only one element different from zero, equal 
to ±1 , in each row and column. The symmetry group is in this case S q x iTLa) q , 
the direct product of S q , the symmetric group of q elements, and (Tto) q ■ 

In the general case (r ^ 0, q), there is no symmetry mixing the black and red 
sectors of the theory. Therefore the symmetry group is simply given by the direct 
product 

G = 0(r) x S t x (Za)* , 

where t = q — r is the number of black indices. 

We now turn our attention to some properties of the transfer matrix r{6) 
which, for a iV -sites system, is defined as 

t(9)=K(R a1 (9)R A2 (9)...R an (9)) . 

From property (2.3) of quasiperiodicity one gets for the transfer matrix 

t(9 + in) = Tr (u A R A1 (9)u A 1 h 1 u A R A 2(9)u A 1 h 2 . . . u A R AN {9)u A A h N ) 
=t(6)H , 

with H = (£)jLi hj . Next we note that the matrix h is a symmetry of R(9) 
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(property d) and that (hu) 1 = u 1 h then 

R 12 (6 + in) = (hu® h)Ri 2 {6){u- l h® 1) . 

So that 

t(9 + in) = Tr [M^ii^X^ 1 ^ {hu) A h 2 R A2 {6){ U - l h) A . . . 
=Ht{6) . 

Thus, 

r 2 (^ + in) = T{e)n 2 r{6) = t 2 {6) (2.4) 

since V? = 1 . This shows that configurations on the square lattice with periodic 
boundary conditions of the NIS model are invariant under the change — > + iit . 

The key property for the applicability of the analytic BA technique to the 
general NIS model is the inversion relation satisfied by the transfer matrix, namely 

r(0 + 77) t{9) = [sinh(?7 + 6) sinh^ - 6)} N + sinh^ 6 t{6) , (2.5) 

where t(8) is a matrix regular at 6 = 0. In order to prove (2.5) it is convenient 
to switch to the S -matrix, which is related to the R -matrix by S^(9) = R^,(d) ■ 
We notice that, from eq. (2.1) and crossing invariance we have 

S$(r]) = smh V 5 a b 6 c d ; 

so S(r]) is proportional to 



1 q 

v = -J2 \ aa ) ( bb \ 



a, 6=1 

which is a projector onto a one-dimensional subspace. The Yang-Baxter equation 
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for the S -matrix reads 



Si 2 (0' - O)S 13 (0)S 23 (0') = S 23 (0')S 13 (0)S 12 (0' - 9) . (2.6) 
Setting 9' — 9 = 7] we get 

V 12 S 13 (9)S 23 (r] + 9) = S 23 ( V + 9)S 13 (9)V 12 
which, left multiplied by = 1 — V\ 2 , yields 

V 12 S 13 (9)S 23 (r ] + 9)Vt 2 = ■ 
Therefore, the matrix S\ 3 {9)S 23 {r] + 9) takes the block triangular form 

S 13 (0)S 23 (r 1 + 9)=[ n ,_J (2-7) 

V * £3(0) / 

in the basis where 7 3 i2 has the canonical form 

*»=(;;)■ ^> 

The matrix notation in eqs. (2.7) and (2.8) refers to the tensor product Vi x V 2 
of the vector spaces V\ and V 2 so that the matrix elements are matrices in the 
vector space V3 . From eq. (2.7) we can determine A{6) which is given by 

V 12 S l3 (9)S 23 ( V + 9)Vi2 = A 3 (9)V i2 ■ 

The ^-matrix defined by the R -matrix elements in eqs. (1.2) is time-reversal in- 
variant 

S(6) = S(9f . 

By using time reversal, unitarity [eq.(2.2)], crossing invariance [ Si 2 (9) T2 = Si 2 {r] — 
9) ] of S(9) and from the explict form of V we get 

A(9)=p(9)l . 

This shows that the matrix A(9) is diagonal. 
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Analogously, we can show that the matrix B{9) is zero at 9 = 0. From 
eq.(2.7) we have 

^^(O)^)^ = £3(0)^2 . (2.9) 
As consequence of crossing and regularity, it is 

where P13 is the exchange operator in the spaces 1 and 3: 

P 13 = \ab) (ba\ . 

a,b—l 

It is then a matter of straightforward algebra to show that the LHS, and therefore 
-B(O) , in (2.9) is zero. With no loss of generality, we can then rewrite eq. (2.7) in 
the form 

( p[Q) \ 

Sl3(0)S23fa + 0)= ■ , fl5 ^ > ( 2 - 10 ) 

V * smh# 03(0 ) J 

where 83(6) is a (q 2 — 1) x (g 2 — 1) dimensional matrix in Vi x V2 regular at 
9 = whose elements are matrices in V3 . 

Consequently 

r(0)r(77 + 6) = Tr (S M (0)S B i{ri + 0)^42(^2(77 + 0) • • • S AN (6)S BN (rj + 0)) 
= p(9) N + sinh^ Tr (£i(0) . . . ^(0)) , 

.423 V / 

with p(0) given in eq. (2.2). This proves eq. (2.5). 
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3. Analytic Bethe Ansatz (BA) for the general NTS model 



We start by constructing a family of eigenstates for the transfer matrix that 
we shall call reference states. Although these states provide only a small portion of 
the spectrum of r{6) , their knowledge is an essential ingredient for the formulation 
of the analytic BA which, in turn, will generate the whole spectrum of t{9) . 

The action of r{6) on a general vector \a\, . . . , a/v) , is given by 
T(9)\ ai ,...,a N )= £ R^mtt::(0)...RZT(0)\b 1 ,...,b N ) ■ (3.1) 

w,m=i 

If we restrict our attention to those states such that [1] a- l+ i ^ a; for i = 
1,...,N, then eq.(3.1) takes a particularly simple form. Indeed, in (3.1) let us 
consider the following two possibilities: (i) c\ ^ a\ and (ii) c\ ^ . In the case 
(i), Rll^(9) is non-zero only if b\ = c\ and C2 = a\ . This implies that in the 
next factor R^ 2 ^(9) , since a\ ^ <i2 , we must have 62 = oi and C3 = a2 and so 
on, hence the sum in (3.1) reduces to just one term: 

N 

]\l%%l\(6)\a N ,a 1 ,a2,...,a N - 1 ) . 
i=i 

Here and in the following we set ajv+i = aj . 

If instead c\ 7^ ajv , then cjv = a-N and c\ = in (9) . In a way 

similar to case (i), this time moving backwards in the product of weight factors in 
eq.(3.1), we get 

N 

n^Wi(^)l a 2---°iV,ai) . 
i=i 

As we have just seen, case (ii) implies a\ — c\ . Thus, the two cases considered 
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above exhaust all the possibilities. The total result is then 



N 



r{9) \ ai ,...,a N ) =\\R a a %l+ + \(6) \a N ,a h a 2 ,...,a N - 1 ) 
i=i 

N 

+ \\K\ + ,a l+1 { e )\ a ^---^N,a l ) . 



i=l 

Now, the value of the weights R^{9) = R^(r] — 9) depends on whether a and b 
are black or red indices. From eqs.(1.2) we find 

Rf h {6) = sinhfo - 9) e* n{a > b) , 

where n(a, b) is the number of black indices among a and b ( n(a, b) =0,1 or 
2). Thus, by taking into account that each index appears twice in the products, 
we obtain 

N N 

U R TaZl( d ) = ^ hN ^-0)e 9B , n^iW) =*nb N eeW 
i=l i=l 

where B is the number of black indices in the state \ai, . . . , ajv) ■ Since the effect 
of applying r{6) amounts, apart from an overall factor, to shifting the indices 
to the left or to the the right, eigenstates of t{9) can be obtained by constructing 
"plane waves" of the form 



L-i „ a. „ 

^(e{;a 1 ,...,a N ) = 2_^e{\a 1+k ,...,a N+k ), (3.3) 

fc=0 a N+i — Ct{ 



with e L = e~ . In eq.(3.3) L is the smallest integer such that . . . , cll+n) — 

\a\, . . . , a/v) and 1 < j < L . Depending on the choice of the indices a\, . . . a/v , 
L takes the values 2 < L < N , with N/L integer. These states are eigenvectors 



12 



of t{6) with eigenvalue 



A(ej, 9) = e{ smh N (r] - 6) e 6B + e~J sinh^ 9 e^ B (3.4) 

and they provide a basis in the sector of the Hilbert space spanned by vectors 
\a\, . . . , ajsf) with Oj+i 7^ Oj . Indeed 

1 L 

|ai, . . . , a N ) = - ai, . . . , a N ) 

owing to Ylf=i £ "l — . Notice that the eigenvalues (3.4) are highly degenerate 
since they do not depend on the particular choice of \a\, . . . , a/v) as long as B is 
left unchanged. 

Clearly B = in the separable case. In the case r = , it is still possible to 
build reference states for which the eigenvalues are of the form (3.4) with B ^ N 
by forming suitable linear combinations of states with one or more couples of equal 
indices. For example starting from states of the form 

\ppp 3 ...p N ) , (3.5) 

with pi 7^ pi + i and p 7^ ps, pjy , as before we form the plane waves 

N-l 

\ £J n^) =^2 £3 N\pi+k,---,PN+k) , (3.6) 

Pi 7^ Pi+l if «7^1, Pl= P2 = P , PN+i = Pi ■ 

Then, states of the kind 

V(e J N ; p, v) = \e J N ; p) - \e J N ; v) 

are eigenstates of the transfer matrix with eigenvalues of the form (3.4) with B = 
N — 2 . By taking appropriate linear combinations of plane waves involving states 
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with two and more couples of equal neighboring indices, we span all the possible 
even values of B from N to (N is supposed to be even). This also shows 
that, in the case where only black indices are allowed, B is bound to take on even 
values. 

We are now in the position to formulate the analytic BA for the general NIS 
model. Eqs. (2.5), (2.4) and crossing invariance of r(6>) imply the following 
properties for the eigenvalues A{9) of r{9) : 

A(0 + 77) A(0) = sinh iV (77 + 9) sinh^ - 9) + sinh^ 9 A{9) , (3.7) 
A(9 + m) 2 = A{9) 2 , (3.8) 



A(9) = A( v -9) , (3.9) 

where A(9) is a regular function at 9 = 0. 

Starting from expression (3.4) of the reference eigenvalues, we propose for the 
A(0) the general analytical form 

A(0) = F{9) sinh^ - 9) e 9B + G{9) sinh^ 9 e^ B , (3.10) 

where B is an integer (0 < B < N ) and F(9) and G(9) are functions to be 
determined. As previously discussed, B = in the separable case and B even 
if r = 0. From the explicit expression of the R -matrix elements and eq.(3.8) 
it follows that A{9) must be a polynomial in e° of degree at most 2A and of 
definite parity under the change 9 — > 9 + iir . Therefore with no loss of generality 
we can take F{9) and G(9) to be rational functions of e 20 with the same number 
of zeros and poles. Moreover, in order to obtain A(0) entire in e° , the poles of 
F(9) and G(9) must have the same location, so that the corresponding singular 
contributions may be made to cancel. This last requirement is what will give us 
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the Bethe-Ansatz equations. Summarizing we look for F{9) and G(9) of the form 



p 26 - 7- 

^=^11^— T' (3.11)(3.11) 

i=i 6 « 



^)=^IIVT' (3.12)(3.12) 

i=i 6 « 



where (f),ip,ZiXi and £j are complex parameters to be determined. The expression 
(3.10) and the form of F(9) and G(9) shows that the parity of A(9) under the 
change 9 — > 9 + in is entirely determined by the choice of B . 

From eqs.(3.7), (3.8) and (3.9) it follows 



F(9) G(9 - rj) = 1, G(9) = F(rj-9) (3.13) 



or 



F(9) F{-6) = 1 . (3.14) 
By substituting the expression (3.11) into eq. (3.14) we obtain 



3 2R e(? !> TT ' ' TT ! ' ■- 



i=l i=l 



Zeros and poles on both sides of eq.(3.15) coincide if ^%Zi = 1. Equating the 
residues at the poles of both members of eq. (3.15) imposes 



m 



On the other hand eqs. (3.13) give 

e ~ ^ —<b TT e 



e 



-i ' 

i=l ~ i=l " s * 



that is: ip = —<p , Q = e 2?? ^ and £j = e 2r? £ • 1 . It is then natural to define 
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Hi — e so that 

Ci = e 3v /j,i, Re(0) = mq with {/^} = j^J 1 ] 
The final result is 



^)=^II^3^' (3.16)(3.16) 

. JUL e 20~V _ p2r?... 

G{0) ^jl ^J^ , (3.17)(3.17) 

with a; = Im0 . 

Reinserting eqs.(3.17) and (3.16) back in (3.10), and imposing the vanishing 
of the residues at the poles of A(6) , so to actually obtain a polynomial in e° , we 
get the Bethe-Ansatz (BA) equations. The singularities of F{9) and G{9) are at 
those values 6k of 9 such that e 26k = e v fik ■ The BA equations are thus: 

sinh 6 k \ N ( v -29 k )B _ 2iuj TT sinh(# fc - ^ + 77) 



n 



sinh(?7 - fc ) / A± s inh(0 fc - ^ - 77) 

which, after the change of variable 6>fc — rj/2 — i\k , read (assuming A even) 

/ sin(A fc + ^/2) \ ^ e 2i(A fc B- w ) = TT sin(A fc - Aj + ig) ^ lg , 

Vsin(A fc - ?V 2 )/ | j- sin(A fe - A; - 177) 

The eigenvalues of the transfer matrix have thus the form 



\J =1 sm(A,- ~ 20 + V2) j w ; 

/" sin(A - l 9 + 3 i r ] /2) \ N ^_ e)B _ lu] 

\H ^(Aj- 20 + 277/2) J 



(3.19) 



— * 

with A = (Ai, . . . , A m ) solution of (3.18). We see that the BA equations of the 
general NIS model are very similar but not identical to those of the six-vertex 
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model. First, in the six-vertex with periodic boundary conditions the phase factor 
e tuJ is absent; second and more important, a A -dependent phase factor e 2l ^ kB 
never appears in the six-vertex model. It is usually possible to determine the 
phase uj by studying the asymptotic behaviour for 9 — > oo of the transfer matrix 
[12]. For the general NIS model the diagonalization of the transfer matrix in this 
limit is not simpler than the diagonalization for finite values of 6 . An analysis of 
the large 9 asymptotics of r(6) will be the subject of the the next section. In the 
separable case e luJ is known to be a root of unity [1]. Computations on small size 
examples and the analysis of the next section led us to conjecture that the same 
result holds true in the general case; in Sec. 6 we provide a proof of this result in 
the thermodynamic limit. 

To analyze the structure of the BA equations, we take as usual the logarithm 
of eq. (3.18) which reads (7 = — 77 > ) 

m 

iV$(A fc , 7/2) + 2X k B -2uj-J2 H^k ~ \, 7) = 27r/ fc (3.20) 

i=l 



where 

sin(ry + A) 



$(A, 7 ) = ilog 



sin(ry — A) 



(3.21) 



and the I k are half integers. We choose the branches of the logarithm in (3.21) in 
such a way that the function $(^,7) for real x is continuous, odd and monoton- 
ically increasing for —it/2 < A < 7r/2 . In particular 



(±|, 7 ) =±vr , (3.22) 



$(A + 7T, 7 ) -$(A,7) = 2tt . (3.23) 
As customary we introduce the so called counting function Zn(\) [11] that in 
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the present case takes the form 



Z N {\) = 



2tt 



$(A, 7 /2) + 2 (bX - - - *(A - A,, 7) 



(3.24) 



where b = B/N is in the range < 6 < 1 . By construction, at each real real root 
Afc Z^(X)N takes the half-integer value (cfr. eqs. (3.20) and (3.24)) 1^ : 



^iv(Afc) 



4 

TV 



(3.25) 



Let us solve the BA equations (3.20) in the thermodynamic limit. For the 
ground state we assume as usual the regular filling 



'fe+i - Jfc 



4 = 1. 



Then, after introducing the continuous density of roots 



a b (X k ) = lim 



iv^oo iV(A fc+ i - A fc ) ' 

in the N — > oo limit in which the BA roots are closely spaced, we have from eq. 
(3.25) [13] 

dZ oo 



dX 



Thus in the limit iV — > oo eq. (3.24) reads 



tt/2 

( 7 6 (A) = i-$ / (A, 7 /2) + ^-^ | ^<t>'(A-//, 7 W//) , 

-tt/2 



(3.26) 



with $'(A, 7) = d$(A, "f)/dX . Eq. (3.26) is identical to the corresponding equa- 
tion for the six- vertex model, except for the constant term b/n [13]. This linear 
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equation can be solved by Fourier series. Starting with 

+00 

a b (X) = ****** 

k——oo 



$(A,7) = 2A + 2^^^e- 2 ^ 
k=l 



we find 

1 + 5 1 +2^ e 2ik\ 

fe^O 

This density results to be equal to the six-vertex ground state density of roots plus 
the constant b/ (2ir) . This implies that the number of roots 

7r/2 

TV ' N + B 



J d\a b (X) 



-7T/2 

is B/2 larger than in the six-vertex model. 

The free energy per site at fixed b is defined as 

A(0, 7 ) = - lim ^logA(tf) 

N^oo JV 

which, by using eq. (3.19), takes the form 

tt/2 

f b (6,>y) = -log|sinh( 7 + 0)| -9b + iJ d\$(\- 16,7/2) a b (\) . (3.28) 

-tt/2 

Inserting eq. (3.27) in (3.28) yields the b -independent result 

Men) = Me,-y) 

in the thermodynamic limit. This shows that all the ground states of the various 
B -dependent sectors of the Hilbert space are degenerate in this limit. 



19 



The excited states follow by allowing holes in the half-integers Ik ■ 

N h 

h+l - h — 1 + ^ Siih ■ 

h=\ 

The corresponding root density and eigenvalues are identical to those of the six- 
vertex model since they are B -independent. 

4. The large 6 behaviour of the NIS 
model and of its Bethe Ansatz solution. 

It is always instructive to investigate the properties of integrable models in 
the limit in which the spectral parameter 9 takes large values. In particular the 
quantum group structure usually emerges in such a limit [13]. We see from eq.(1.2) 
that the leading contribution for large 9 of the various Bolztmann weights is e de 
with d = 2, |,l,i,0. This behaviour makes the diagonahzation of the transfer 
matrix in the large 9 limit as difficult as for finite values of 9 . This situation 
differs from the one of the six- vertex model where d — 1, and the transfer matrix 
in the limit 9 — > oo is easily diagonalizable and non-singular. 

Clearly, the purely black configurations R^u(9) = R^(r] — 9) [see eqs.(1.2)] 
dominate in the 9 — > +oo limit. More precisely 

RZ(°) „ = e 2 VMi - <W) • 

>+oo z 

One then gets for the N -sites transfer matrix in the same limit 

/ -v\ N 

'<v^(-V) +*-♦>) ' 

where 

M 2N = ^ I//1//2 ■ ■ ■ m) (^2/^3 ■ ■ ■ A*iVA*l| , Vi^Vi+i (4-1) 

is a left-shift operator in the purely black sector with the condition that all the 
neighboring indices are different. Notice that the operator M 2N has a non-zero 
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kernel since it annihilates all states having either a couple of equal neighboring 
black indices or at least one red index. Thus the diagonalization of M 2N does not 
yield any information on the asymptotic behaviour at large 9 of the states in its 
kernel. We then move to the next to leading contribution in r(9) . It is obtained 
by replacing in two neighboring sites the weights R l j Z(9) with R^a(9) and Ra^,{9) 
in all possible ways. From 



Km = Km 



e2 c 



we find 



N 

(n o ^ —J (e 2m M 2N + e^ N -^M 2N _) (l + o(e^)) , (4.2) 



where the operator M 2N _ 1 is again a left-shift operator of the form (4.1) with all 
black indices but one which is red. Just like M 2N , the operator M 2N _ 1 has a 
non-zero kernel. 

It is useful to match eq. (4.2) with the large 9 behaviour of the eigenvalues 
(3.19) obtained from the BA equations. We have 



A(9) 



++oo \ 2 



N 



M JN+B)0 .iu-(N-m)v 



l + e- 26 2sinh??^e 



-2i\j _ j\f e 2n 



3=1 



N 



_|_ I _ ] e (N-B)9 e -iLU+(B-m) V 



+ lower orders in e 



1 - e- 29 [ 2e 2v sinh r] e" 2iAj + N 

3=1 

(4.3) 

Comparison between eqs.(4.3) and (4.2) leads to the conclusion that the highest 
degree term e 2Ne in eq. (4.2) corresponds to an eigenvector ip 2N of M 2N whose 
eigenvalue in the large 9 limit is given by (4.3) with B = N . The eigenvectors of 
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M 2 „ are reference states of the form (3.3) with all black indices. Moreover 

M 2N ^j 2N = e{ip N 

since M 2N is a shift operator. Thus from eq. (4.3) e tuJ must be a L-th root of 
unity. The operator M 2N annihilates all vectors outside the purely black sector 
of the Hilbert space. States with only one red index are instead captured by the 
operator M 2N _ 1 which, being a shift operator, has eigenstates of the form (3.3) 
and roots of unity as eigenvalues. Thus also in this case e lu} is a root of unity. 
Again there are no BA roots if B — N — 1 . 

In order to find roots of the BA equations we must move to the next to sub- 
leading term, e^ 2N ~ 2 ^ , corresponding to B = N — 2 in the asymptotic expansion 
of t{6) for large 9 . The transfer matrix in this limit is 

/ -v\ N 

r{0) = (-V) \e 2N6 M 2N + e^ N -^M 2N _ 1 + e^ N -^M 2N J 



+ 



—M 2N + M 



The operator M 2N _ 2 is a shift operator of the form (4.1) involving 2 red indices. 
Its eigenstates are then of the kind (3.3) and its eigenvalues are roots of unity. The 
operator M is given by 

M = sinh V M^ 2 + yJl41 + ^M 2 ( i 2 + 1 ~ 6 2 — Mff_ 2 + l -M 2 ^}_ 2 + ^M 2 ( ?_ 2 

where the operators M 2 l }_ 2 ( i = 1, . . . , 6 ), as the diagrammatic representation of 
in Fig. 2 shows, involve a couple of equal neighboring indices. Eigenvalues of the 
transfer matrix with asymptotic behaviour e( 2N ~ 2 ^® not corresponding to reference 
states can then be obtained by diagonalizing the operator M . Obviously (see 
Fig. 2) the reference states (3.3) belong to the kernel of M and conversely the 
eigenstates of M are in the kernel of M 2N , M 2N ^ and M 2JV _ 2 . Let us start by 
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considering the plane waves 



e J N ; nj of eq.(3.6). Calling and vr the (black) 

indices on the left and on the right of the couple /ifj, in the states sP N \ (xj , the 
state 



*(4) = e" /2 H^) + 2e- r '/ 2 smhr ] \ei-,u R ) + J2H-,a) 

m=i a=l 



satisfies 



where eq. (1.3) has been taken into account. Therefore we have 

This corresponds to a solution of the BA equations with B = N — 2 , m = 2 and 
e iW = £^r which is a root of unity. 

In conclusion, we have verified for B = N, N — 1, N — 2 that e laj is a root of 
unity. We conjecture that this is the case for all values of B . 

5. The partition function and the excitation 
spectrum in the thermodynamic limit 

In this section we collect some results about the determination of the free 
energy per site in the thermodynamic limit and make the comparison between the 
result obtained by means of inversion techniques [2,14-16] and the one we get from 
the BA equations (3.18). 

We start by observing that eq.(2.5) can be rewritten as 

r(e + r ] )T(e)=p(6) N l + o(e- N ) , (5.1) 

where now we have made explicit the fact that the second term on the RHS is 
vanishing in the limit iV — > oo . Correspondingly, for the eigenvalues A(9) of t{6) 
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one gets 



lim A(6)A(ri + 6) = p(6) N . (5.2) 

AT— >oo 



For the ground state k${6) we define 

H&) = I™ Ao(0) 1/N 

sinn 77 N^oo 

which, because of (3.9) and (5.2), satisfies (Aq(6>) is real) 



sinh(r] + 6) sinh(^ — 6) 
sinh 2 77 



l q {0)H-o) = - vVM-e) , (5.3) 



where property (3.8) has been taken into account and a suitable function tp(0) 
has been introduced. The explicit form of ip{6) has to be chosen in such a way 
that the following properties for Lq{9) are satisfied: i) analyticity and absence of 
zeroes in the physical strip — r\ < Re9 < 77 , ii) in periodicity, which follows from 
eq. (3.8) under the hypothesis that the ground state belongs to the even B sector 
as discussed in the Sec. 3. 

Eq.(5.3) can be solved iteratively. Starting from Lq{9) = <p(9) and imposing 
alternatively crossing invariance (3.9) and eq.(5.3) we end up with the solution in 
terms of an infinite product 

L (g) = TT gffi?? + 9) tPi^k + Vv-e) <p(2(k + l) V ) 
0U ll Q <p((2k + l)7) + 9) <p(2(k + l)7)-9) V {2k v ) 

where the last ratio in the product has been introduced so as to make Io(0) = 
I/O (77) = 1 . The obvious identification is then 

sinhfa + 9) 

Pi u = 

smhTj 

so that, in the thermodinamic limit, 

\(0\V N = --h TT sinh (( 2fc + ^ + g) sinh(2(A; + l)rj - 9) sinh((2A: + 3)77) 
01 ' " Sm V 11 sinh(2(A; + 1)77 + 9) sinh((2A; + 3)t? - 6) sinh((2A; + 1)77) ' 

(5.4) 

This solution has the required analyticity and periodicity properties. By taking 
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the logarithm of Aq(9) we obtain the free energy per site fms(9,v) 



f ms (9,r ] ) = -\ogA (9) 1 / N . 

The expression (5.4) for the ground eigenvalue could have been obtained equally 
well from the BA equations (3.18) with B = . From eq.(3.28) we get for the free 
energy 

/m 4»,, ) ) = » + EV^W- log i stah (''- (,)l ' (5 ' 5) 

k — 1 

The sum on the RHS can be recast in the following form 

y> sinh(2fcfl) _ y> fc 1 - e 2(g+fc^) 
cosher?) ~ ^ ( ^ ° g i _ e -2(9-k v ) 

k — 1 k- — 1 

1 _ e -29+4kri ]_ _ e 2^+2(2fc-l)77 



-logllx 



e -26»+2(2fc-l)r? X - e 2e+Akr > 

k=\ 



so that (see also [16]) 

A Q (9) l ' N =e-Sm 

00 1 _ e -26+4(k+l) V l _ e 29+2{2k+l)rj 



2 Wi_ e -s 



2 11 1 _ e -26»+2(2fc+3)r; ]_ _ e 26>+4(fc+l)?7 
fc=0 



which is completely equivalent to (5.4) owing to 



oo 



^-r sinh((2fc + 3)7y) = 

sinh((2/c + l)rj) 2sinh?7 

k — 

The derivation above shows the relation between the free energies of the NIS and 
six vertex models is 

/nis(M) = /e-v(-0,7) » 1 = -V ■ 

Following ref. [15], the inversion method can also be used for determining all 



25 



those eigenvalues A(9) such that the limit 



w = £!Lm (5 ' 6) 

is finite and non-zero. From the definition (5.6) and eq.(5.2), 1(9) satisfies the 
inversion relation 



l(9)l(rj + 9) = 1 or l(9)l(-9) = 1 , (5.7) 

where in the last equation crossing invariance has been taken into account. From 
crossing invariance and eq. (5.7) it is immediate to show that 1(9) , in addition to 
being in periodic owing to property (3.8), is also 2r\ periodic: 

1 



1(9) = l( v -9) = — = = = 1(9 - 2rj) . 

1(0 -v) l(2r]-9) 

The excitations, just like the ground state, ought to be analitic in the physical strip 
— rj <Re9<r] but zeros are now allowed. The typical ansatz is thus written as a 
product of Jacobi elliptic functions of modulus k as follows 

n 

l(9) = l[C J sn(D(9 J - i 9)) , 

with Cj and D coefficients to be determined. In the following we shall assume n 
to be even, as in the six-vertex model. The parameters 9j specify the location of 
the zeroes. Periodicity 2i] and m require 

K' , K' 2ri , . 

D — — and — = . 5.8 

1] K n 



By making use of [17] 



sn(^ — iK') sniz) = — , (5.9) 
k 



2G 



eqs.(5.7) are satisfied if 



Cj = ±Vk , and {Oj} = {0i - irj} 



The simplest solution, corresponding to elementary excitations, is given by 



9j = otj + 277/2 , f < j < n 



(5.10) 



where real numbers (other more general solutions solutions would describe 
bound states). For the given by eq.(5.10) we get 



n f K' K'\ 

1(9) = ± Y[ Vk sn — (aij - i&) + i— 



(5.11) 



The energy per site of the ground state is given by 



eo = 



df NIS (0) 



d9 



2 E 

fc=i 



cosh(krj) sinhr/ 



The momenta and energies of the first excitations can be obtained from the ex- 
pression of 1(9) . We have 



P-Po = -ilnZ(O) = ^Pj , 
d\nl(9) 



e - e = - 



d9 



■3} ' 



=0 , =1 



(5.12) 



which give 



,K' K' 
Pj = —tvk logsn — a* +t — 
v 77 2 



(5.13) 
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and, with the help of [17], 



d snu 




2 u){l — k 2 sn 2 u) , 



du 



one gets 




(5.14) 



Notice that pj is real owing to eq. (5.9) as it must be. From eq. (5.14) the energy 
gap is obtained by setting pj = with the minimal value n = 2 in eqs.(5.12): 



In summary, the free energy (5.5) and the excited states (5.11) of the general 
NIS model in the thermodynamic limit are completely equivalent to those of the 
six-vertex ones. There is one branch of excitations with non-zero gap, given by eq. 
(5.15), and obeying the dispersion relation (5.14). 

Finite size corrections can be computed here with the method proposed in [18]. 
Since the model is gapful, they are exponential in the size. More precisely they are 
typically of order 




(5.15) 



k+ N / 2 , 



where k < 1 is defined by eq. (5.8). 
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6. The NIS model as a mixed spin 1/2 and oo-spin model 



In this section we show that the BA equations and the transfer matrix eigen- 
values of the NIS model derived in Section 3 can be identified with those of a spin 
system describing a mixture of spin 1/2 and spin s particles in the limit in which 
s — > oo . In ref. [10] it has been shown that the ground state of a pure spin s 
model is formed by a Dirac sea of solutions of the BA equation arranged in strings 
of 2s roots (called 2s -strings). This situation generalizes the usual situation en- 
countered in spin 1/2 systems where the sea is formed by real roots (1-strings). 
It is also known [9] that the ground state of a system with mixed spin 1/2 and 
s — 1 is just the superposition of the two seas. Since the form of the BA equations 
found in [9] is immediately generalizable to an arbitrary value s , we can conclude 
that the ground state of a system with mixed 1/2 and s spin is formed by a sea of 
2s -strings plus a sea of real roots. Namely, we have M real roots Xj ( 1 < j < M ) 
and v s 2s -strings: 



K,l = cr a + i(l + \)^ + e Nl/2Ns 



with 1 < a < u s and — s < I < s — 1 . Here the real numbers cr a are the strings' 
centers and £n 1/2 n s are quantities vanishingly small when N s — > oo . 

The BA equations (with periodic boundary conditions) in presence of string solu- 
tions are then 



/ sin(A fc + ? 7 /2) \ Nl ' 2 / sin(A fc + ?s 7 ) \ ^ = _ TT sin(A fc - A, + 27) 
Vsin(A fc - ry/2)/ \sin(\ k - is^) J Al sm(X k - A; - ij) 

x nn sin(A fc - X aJ + ij) 
a=ij=- a sin(A fc - A aJ - i 7 ) 



(6.1) 
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for a real root A& and 



sin 



Vsin(A Q)/ - 17/2)/ \ sin ( A a,; - is7)/ sin(A ai ; - A; - 17) 

x TT TT sin(A Q; ; - A Aj + yy) 
s sin (A Q ,Z - A/Jj - 17) 
(6.2) 

for a string solution A a ^ . The integers jV]y 2 and N s denote the number of sites 
where particles with spin 1/2 and s respectively sit. Both N1/2 and N s are 
assumed to be even integers. 

The transfer matrix eigenvalue of this mixed spin system is given by 

A(u, A) = Aa(u, A) + A-£,(u, A) , (6.3) 

with 

A A (u, A) = sinh^ (7 _ u) ( smh[(, + 1/2)7 - ,] \ * rr sin^ + iu + » 7 /2) 
V ' ; u ; V sinh[(s + 1/2)7] / ^ sin(A j + iu - i 7 /2) 

rj iht sin(A Q; ; + iu + i^/2) 

(6.4) 

— * 

where A = (Ai, . . . , A^, Ai,- S , . . . , X Vs , s -i) are solutions of (6.1) and (6.2) and 

— * ~ 

A D (u, A) = A A (~f-u, A). 

The infinite s limit of eq. (6.1) gives then 

( sm(A fc + 2 7 /2) \ Nl " -2i{N.-2u.)\ k _ _ pxn V ff ^ TT sin ( A fe~ A » + ^) 
Uin(A,-i 7 /2)J " exp ( v 4 ^ i (T ^ll s in(A fe -A i -i 7 ) 

(6.5) 

which coincide with eqs.(3.18) if we set 7 = — i] , = AT, B = N s — 2v s and 

e «o; _ eX p(— 2i X)a=i °a) ■ ^ n ^ ne same wa Y an( i with the same identification of 
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parameters the large s limit of eq.(6.4) is 



A A (u,\) = sinh iVl / 2 (7 - u)e 



a=l 



-Mj sin(Aj + ra — ry/2) 



which inserted back in (6.3) reproduces exactly the eigenvalue (3.19) of the NIS 
model if u = —9 . 

One can further explore this correspondence and give a proof, at least in the 
limit iV — > oo , to our conjecture that the phase factor e tuJ of the NIS model is 
in fact a root of unity. To this purpose we multiply all the BA equations (6.2) 
corresponding to a string with center say a a and take the logarithm [19]. After 
some algebra we obtain in the limit Ny2,N s — > oo 



2s 



fc=l /3=1 



$(a Q - ap, 2s 7 ) 



2s-l 



+2 Qioa-Ofaki) 
fc=l 

where, for real A , 



+ E [ $ " ( s + sb) + $ (<t« - A„ (s - i) 7)] 

(6.6) 



<3>(A, 7 ) = 2tan 1 [tan A coth 7 ] . 

As usual I a denotes a half-integer. Eqs. (6.6) are a set of equations for the real 
roots and the centers of the strings. We can easily solve for the latter in the infinite 
s limit in which (6.6) reduces to 

v s M 
{N 1/2 + 4sN s )a a = 2nl a + 8s ^(a Q - ap) + 2 ^(a a - \j) + o(s°) . 

(3=1 j=\ 

Under the hypothesis of regular filling for the ground state 

N- M 



M=^ + o(l), J> = (1) 
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we get 




This result shows that the quantity 



cxp 



2 *£ 

0=1 



Or 



cxp 



ITX 

7n, 



/3=i 



appearing in (6.5) is indeed a root of unity and the correspondence with the eigen- 
values of the NIS model holds. 

Recall that the eigenvalue of S z on a state determined by eqs. (6.1)-(6.4) is 
given by 

S z = s(N s - 2u 8 ) + - M 

The first term is the contribution of the s spins and the second one of the spins 
1/2. Then if we define 

B — lim — . 

s— »oo s 

Hence B can be interpreted as a properly normalized eigenvalue of S z for s = oo . 
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FIGURE CAPTIONS 

1) Boltzmann weights for the general NIS model. Dashed and solid lines repre- 
sent red and black indices. The remaining non-zero weights can be obtained 
by crossing invariance. 

2) Contribution of order e ( 2N ~ 2 )° to the transfer matrix in the limit 6 — > oo . 
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